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ABSTRACT 

Sun-like stars have stellar, brown dwarf and planetary companions. To help constrain their formation 
and migration scenarios, we analyse the close companions (orbital period < 5 years) of nearby Sun- like 
stars. By using the same sample to extract the relative numbers of stellar, brown dwarf and planetary 
companions, we verify the existence of a very dry brown dwarf desert and describe it quantitatively. 
With decreasing mass, the companion mass function drops by almost two orders of magnitude from 
1 Mq stellar companions to the brown dwarf desert and then rises by more than an order of magnitude 
from brown dwarfs to Jupiter-mass planets. The slopes of the planetary and stellar companion mass 
functions are of opposite sign and are incompatible at the 3 sigma level, thus yielding a brown dwarf 
desert. The minimum number of companions per unit interval in log mass (the driest part of the 
desert) is at M = 31 j^^g Mjup- Approximately 16% of Sun-like stars have close (P < 5 years) 
companions more massive than Jupiter: 11% ±3% are stellar, < 1% are brown dwarf and 5% ±2% are 
giant planets. The steep decline in the number of companions in the brown dwarf regime, compared to 
the initial mass function of individual stars and free-floating brown dwarfs, suggests either a different 
spectrum of gravitational fragmentation in the formation environment or post-formation migratory 
processes disinclined to leave brown dwarfs in close orbits. 

Subject headings: 



1. INTRODUCTION 

The formation of a binary star via molecular cloud 
fragmentation and collapse, and the formation of a mas- 
sive planet via accretion around a core in a protoplan- 
etary disk both involve the production of a binary sys- 
tem, but are usually recognized as distinct processes (e.g. 
Heacox 1999; Kroupa & Bouvier 2003, see however Boss 
2002). The formation of companion brown dwarfs, with 
masses in between the stellar and planetary mass r anges, 
may have elements of both or some new mechanism ijBatd 
11)00; Rice et al. 2003; Jiang. Laughlin & Lin 2004- For 
the purposes of our analysis brown dwarfs can be con- 
veniently defined as bodies massive enough to burn deu- 
terium (M > 13 Mjup), but not massive enough to burn 
hydrogen {M < 80M,/„p e.g. Burrows 1997). Since fusion 
does not turn on in gravitationally collapsing fragments 
of a molecular cloud until the final masses of the frag- 
ments are largely in place, gravitational collapse, frag- 
mentation and accretion should produce a spectrum of 
masses that does not know about these deuterium and 
hydrogen burning boundaries. Thus, these mass bound- 
aries should not necessarily correspond to transitions in 
the mode of formation. The physics of gravitational 
collapse, fragmentation, accretion disk stability and the 
transfer of angular momentum, should be responsible for 
the relative abundances of objects of different masses, 
not fusion onset limits. 

However, there seems to be a brown dwarf desert - a 
deficit in the frequency of brown dwarf companions ei- 
ther relative to the frequency of less massive planetary 



companions l|Marcv &: Butleil2000() or relative to the fre- 
quency of more massive stellar companions to Sun-like 
hosts. The goal of this work is (i) to verify that this 
desert is not a selection effect due to our inablility to de- 
tect brown dwarfs and (ii) to quantify the brown dwarf 
desert more carefully with respect to both stars and plan- 
ets. By selecting a single sample of nearby stars as po- 
tential hosts for all types of companions, we can better 
control selection effects and more accurately determine 
the relative number of companions more and less massive 
than brown dwarfs. 

Various models have been suggested for the formation 
of companion stars, brown dwarfs and planets (e.g. Lar- 
son 2003, Kroupa & Bouvier 2003, Bate 2000, Matzner 
& Levin 2004, Boss 2002, Rice et al. 2003). AU mod- 
els involve gravitational collapse and a mechanism for the 
transfer of energy and angular momentum away from the 
collapsing material. 

Observations of giant planets in close orbits have chal- 
lenged the conventional view in which giant planets form 
beyond the ice zone and stay there (e.g. Udry 2003). 
Various types of migration have been proposed to meet 
this challenge. The most important factors in determin- 
ing the result of the migration is the time of formation 
and mass of the secondary and its relation to the mass 
and time evolution of the disk (e.g. Armitage & Bonnell 
2002). We may be able to constrain the above models 
by quantitative analysis of the brown dwarf desert. For 
example, if two distinct processes are responsible for the 
formation of stellar and planetary secondaries, we would 
expect well-defined slopes of the mass function in these 
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Fig. 1. — Our Close Sample. Hertzsprung-Russell diagram 
for Hipparcos stars closer than 25 pc. Small black dots are 
Hipparcos stars not being monitored for possible companions 
W one of the 8 high precision Doppler surveys considered here 
JLineweaver & Grethcr 2003). Larger blue dots are the subset of 
Hipparcos stars that are being monitored ("Target Stars") but have 
as yet no known planetary companions. The still larger red dots 
are the subset of target stars hosting detected planets ("Planet 
Host Stars") and the green dots are those hosts with larger mass 
(M2 > ISMjup) companions ("Other Host Stars"). Only compan- 
ions in our less-biased sample (P < 5 years and M2 > 10~^Mq) 
are shown (see Section 2.2). Our Sun is shown as the black cross. 
The grey parallelogram is the region of My - (B — V) space that 
contains the highest fraction (as shown by the triangles) of Hippar- 
cos stars that are being monitored for exoplanets. This Sun-like 
region - late F to early K type main sequence stars - contains 
our Hipparcos Sun-like Stars. The target fraction needs to be as 
high as possible to minimize selection effects potentially associated 
with companion frequency. The target fraction is calculated from 
the number of main sequence stars, i.e., the number of stars in 
each bin between the two dashed lines. This plot contains 1509 
Hipparcos stars, of which 627 are Doppler target stars. The Sun- 
like region contains 464 Hipparcos stars, of which 384 are target 
stars. Thus, the target fraction in the Sun-like grey parallelogram 
is ~ 83% (= 384/464). 

mass ranges to meet in a sharp brown dwarf valley. 

We examine the mass, and period distributions for 
companion brown dwarfs and compare them with those 
of companion stars and planets. The work most similar 
to our analysis h as been carried out b y iHeacqx. 11999. ): 
Zucker fc Mazehl ll2001bD and Mazch ^alr |2003ji . 
Heacoxl l|1999l) and lz7cker fc Mazeh (,.200Tbl) both com- 
bined the stellar sample of iDuauennov fc Mavod l)1991|) 
along with the known substellar companions and iden- 
tified different mass functions for the planetary mass 
regime below 10 Mj^p but found similar flat distribu- 
tions in logar ithmic mass fo r brown dwarf and stellar 
companions. iHeacod l|1999() found that the logarith- 
mic mass function in the planetary regime is best fit 
by a power-law with a slightly negative slope whereas 
iZucker fc Maze h (2001 11) f o und a n approximately flat 
distribution. [Mazeh et ol\ l)2003j) looked at a sample 
of main sequence stars using infrared spectroscopy and 
combined them with the known substellar companions 
and found that in log mass, the stellar companions re- 
duce in number towards the brown dwarf mass range. 
They identify a flat distribution for planetary mass com- 
panions. We discuss the comparison of our results to 
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Fig. 2. — Our Far Sample. Same as Fig. 0but for all Hippar- 
cos stars closer than 50 pc. The major reason the target fraction 
(~ 61%, triangles) is lower than in the 25 pc sample (~ 83%) is 
that K stars become too faint to include in many of the high pre- 
cision Dop pler surveys where the apparent magnitude is limited 
to V < 7.5 iLineweaver fc Greth'eri>200aV This plot contains 6924 
Hipparcos stars, of which 2351 are target stars. The grey paral- 
lelogram contains 3296 Hipparcos stars, of which 2001 are high 
precision Doppler target stars (61% ~ 2001/3296). The stars be- 
low the main sequence and the stars to the right of the M dwarfs 
are largely due to uncertainties in the Hipparcos parallax or B — V 
determinations. 

these in Section 3.1. 

2. DEFINING A LESS BIASED SAMPLE OF COMPANIONS 

2.1. Host Sample Selection Effects 

High precision Doppler surveys are monitoring Sun-like 
stars for planetary companions and are necessarily sensi- 
tive enough to detect brown dwarfs and stellar compan- 
ions within the same range of orbital period. However, to 
compare the relative abundances of stellar, brown dwarf 
and planetary companions, we cannot select our poten- 
tial hosts from a non-overlapping union of the FGK spec- 
tral type target stars of the longest running, high preci- 
sion Doppler surveys that are being monitored for planets 
ijLineweaver fc Gretheil l2?)()^ . This is because Doppler 
survey target selection criteria often exclude close bina- 
ries (separation < 2") from the target lists, and are not 
focused on detecting stellar companions. Some stars have 
also been left off the ta rget lists because o f high stellar 
chromospheric activity ijFischer et a/.lll'999j) . These sur- 
veys are biased against finding stellar mass companions. 
We correct for this bias by identifying the excluded tar- 
gets and then including in our sample any stellar compan- 
ions from other Doppler searches found in the literature. 
Our sample selection is illustrated in Fig. ^and detailed 
in Table ^ (complete list in the electronic version only) 
for stars closer than 25 pc and Fig. [3 for stars closer 
than 50 pc. 

Most Doppler survey target stars come from the Hip- 
parcos catalogue because host stars need to be both 
bright and have accurate masses for the Doppler method 
to be useful in determining the companion's mass. One 
could imagine that the Hipparcos catalogue would be 
biased in favor of binarity since hosts with bright close- 
orbiting stellar companions would be over-represented. 
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Fig. 3. — Fraction of stars that are known to be close (P < 5 
years) Doppler binaries as a function of absolute magnitude. For 
the 25 pc Sun-like sample (large dots), ~ 11% of stars are binaries 
and within the error bars, brighter stars do not appear to be sig- 
nificantly over-represented. If we include the extra stars to make 
the 50 pc Sun-like sample (small dots), the stellar binary fraction 
is lower and decreases as the systems get fainter. 
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Fig. 4. — Distance Dependence of Sample and Companions. Here 
we show the number of nearby Sun-like stars as a function of dis- 
tance. Each histogram bin represents the stars in an equal volume 
spherical shell. Hence, a sample that is complete in distance out 
to 50 pc would produce a flat histogram (indicated by the horizon- 
tal dashed line). The lightest shade of grey represents Hipparcos 
Sun-like Stars out to 50 pc that fall within the parallelogram of 
Fig. 121 ("HSS"). The next darker shade of grey represents Hip- 
parcos stars that are being monitored for planets using the high 
precision Doppler techniques (8 groups described in Lineweaver & 
Grether 2003). The triangles represent this number as a fraction 
of Hipparcos stars. This fraction needs to be as large as possible 
to minimize distance dependent selection effects in the target sam- 
ple potentially associated with companion frequency. Also shown 
(darker grey) are the number of Hipparcos stars that have one 
or more companions in the mass range 10""^ < M/Mq < 1, and 
those that host planets (darkest grey). Only those companions in 
the less-biased sample, P < 5 years and M2 > 10~^Mq are shown 
(Section 2.2). The fraction of stars having an apparent magnitude 
V brighter (lower) than a given value are shown by the 5 dotted 
lines for V < 7.5 to y < 9.5. 

We have checked for this over-representation by looking 
at the absolute magnitude dependence of the frequency 
of stellar binarity for systems closer than 25 and 50 pc 



(Fig. 0). We found no significant decrease in the fraction 
of binaries in the dimmer stellar systems for the 25 pc 
sample and only a small decrease in the 50 pc sample. 
Thus, the Hipparcos catalogue provides a good sample 
of potential hosts for our analysis, since it (i) contains 
the Doppler target lists as subsets (ii) is volum e-limited 
for Sun- like stars out to ~ 25 pc ()Reidll200^ and (iii) 
it allows us to identify and correct for stars and stellar 
systems that were excluded. We limit our selection to 
Sun-like stars (0.5 < B — V < 1.0) or approximately 
those with a spectral type between F7 and K3. Follow- 
ing Udry (private communication) and the construction 
of the Coralie target list, we limit our anaylsis to main se- 
quence stars, or those between -0.5 and -1-2.0 dex (below 
and above) an average main sequence value as defined 
by 5.4(B -V) + 2.0 < M„ < 5.4(B - V) - 0.5. This 
sampled region, which we will call our "Sun-like" region 
of the HR diagram, is shown by the grey parallelograms 
in Figs. 

The Hipparcos sample is essentially complete to an ab- 
solute visual magnitude of My = 8.5 (Reid 2002) within 
25 pc of the Sun. Thus the stars in our 25 pc Sun-like 
sample represent a complete, volume-limited sample. In 
our sample we make corrections in companion frequency 
for stars that are not being targeted by Doppler surveys 
as well as corrections for mass and period companion de- 
tection selection effects (see Section 2.2). The result of 
these corrections is our less-biased distribution of com- 
panions to Sun-like stars within 25 pc. We also analyse 
a much larger sample of stars out to 50 pc to understand 
the effect of distance on target selection and companion 
detection. Although less complete, with respect to the 
relative number of companions of different masses, the 
results from the 50 pc sample are similar to the results 
from the 25 pc sample (Section 3). 

Stars in our Sun-like region are plotted as a function 
of distance in Fig. 0] Each histogram bin represents an 
equal volume spherical shell hence a sample complete in 
distance would produce a flat histogram. Also shown are 
the target stars, which are the subset of Hipparcos stars 
that are being monitored for pla nets by one of the 8 high 
precision Doppler surveys (Line weaver fc Gretheill2003j) 
analysed here. The triangles in Fig. 0| represent this 
number as a fraction of Hipparcos stars. 

Since nearly all of the high precision Doppler sur- 
veys have apparent magnitude limited target lists (often 

V < 7.5), we investigate the effect this has on the total 
target fraction as a function of distance. The fraction of 
stars having an apparent magnitude V brighter (lower) 
than a given value are shown by the 5 dotted lines for 

V < 7.5 to V < 9.5. For a survey, magnitude limited to 

V = 7.5, 80% of the Sun-like Hipparcos stars will be ob- 
servable between pc and 25 pc. This rapidly drops to 
only 20% for stars between 48 and 50 pc. Thus the ma- 
jor reason why the target fraction drops with increasing 
distance is that the stars become too faint for the high 
precision Doppler surveys to monitor. The fact that the 
target fraction (triangles) lie near the V < 8.0 line in- 
dicates that on average V ^ 8.0 is the effective limiting 
magnitude of the targets monitored by the 8 combined 
high precision Doppler surveys. 

In Fig. HI 80(= 464 - 384) or 17% of Hipparcos Sun- 
like stellar systems are not present in any of the Doppler 
target lists. The triangles in Fig. ^ indicate that the 
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TABLE 1 
Sun-like 25 pc Sample 



Hipparcos 


B-V 


My 


Distance 


Exoplanot 


Companion 


Number 






(pc) 


Target 


(P < 5 years) 












(M > Mj^p) 


HIP 171 


0.69 


5.33 


12.40 


Yes 




HIP 518 


0.69 


4.44 


20.28 


No 


Star 


HIP 544 


0.75 


5.39 


13.70 


Yes 




HIP 1031 


0.78 


5.68 


20.33 


Yes 




HIP 1292 


0.75 


5.36 


17.62 


Yes 


Planet 



Note. — Table is published in its entirety in the electronic edition of the 
Astrophysical Journal. A portion is shown here for guidance regarding its form 
and content. 



ones left out are spread more or less evenly in B-V space 
spanned by the grey parallelogram. Similarly in Fig. |21 
1295(= 3296 - 2001) or 39% are not included in any 
Doppler target list, but the triangles show that more K 
stars compared to FG stars have not been selected, again 
pointing out that the lower K dwarf stellar brightness is 
the dominant reason for the lower target fraction, not an 
effect strongly biased with respect to one set of compan- 
ions over another. 

In the Sun-like region of Fig. ^we use the target num- 
ber (384) as the mother population for planets and brown 
dwarfs and the Hipparcos number (464) as the mother 
population for stars. To achieve the same normaliza- 
tions for planetary, brown dwarf and stellar companions 
we assume that the fraction of these 384 targets that 
have exoplanet or brown dwarf companions is represen- 
tative of the fraction of the 464 Hipparcos stars that 
have exoplanet or brown dwarf companions. Thus we 
renormalize the planetary and brown dwarf companions 
which have the target sample as their mother population 
to the Hipparcos sample by 464/384 — 1.21 ("renormal- 
ization" ) . Since close-orbitting stellar companions are 
anti-correlated with close-orbitting sub-stellar compan- 
ions and the 384 have been selected to exclude separa- 
tions of < 2", the results from the sample of 384 may 
be a slight over-estimate of the relative frequency of sub- 
stellar companions. However, this over-estimate will be 
less than ^ 11% because this is the frequency of close- 
orbitting stellar secondaries. 

A non-overlap ping sample of the 8 high precision 
Doppler surveys ijLineweaver fc Gretheill2b03|l is used as 
the exoplanet target list where the Elodie target list was 
kindly provided by C. Perrier (private communication) 
and additional information to construct the Coralie tar- 
get list from the Hipparcos catalogue was obtained from 
S. Udry (private com munication). The Ke ck and Lick 
target lists are t hose oflNidever al\ l)2002() . since ~ 7% 
of the targets inl Wright et all ll2004|) have not been ob- 
served over the full 5 year baseline used in this analysis. 
For more details about the sample sizes, observational 
durations, selectio n criteria and sensitivities of th e 8 sur- 
veys see Table 4 of E meweaver fc Gretherill200l . 

2.2. Companion Detection and Selection Effects 

The companions to the above Sun-like sample of host 
stars have primarily been detected using the Doppler 
technique (but not exclusively high precision exoplanet 
Doppler surveys) with some of the stellar pairs also being 



detected as astrometric or visual binaries. Thus we need 
to consider the selection effects of the Doppler method 
in order to define a less-bia sed sample of companions 
(jLineweaver fc Gretheill200^ . As a consequence of the 
exoplanet surveys' limited monitoring duration we only 
select those companions with an orbital period P < 5 
years. To reduce the selection effect due to the Doppler 
sensitivity we also limit our less-biased sample to com- 
panions of mass M2 > O.OOlAf0. 

Fig. [5I shows all of the Doppler companions to the 
Sun-like 25 pc and 50 pc samples within the mass and 
period range considered here. Our less-biased compan- 
ions are enclosed by the thick solid rectangle. Given a 
fixed number of targets, the "Detected" region should 
contain all companions that will be found for this re- 
gion of mass-period space. The "Being Detected" region 
should contain some but not all companions that will be 
found in this region and the "Not Detected" region con- 
tains no companions since the current Doppler surveys 
are either not sensitive enough or have not been observ- 
ing for a long enough duration. To avoid the incomplete 
"Being Detected" region we limit our sample of comp an- 
ions to M2 > O.OOIMq. In iLineweaver et all (2003) we 
describe a crude method for making a completeness cor- 
rection for the lower right corner of the solid rectangle 
falling within the "Being Detected" region. The result 
for the d < 25 pc sample is a one planet correction to 
the lowest mass bin and for the d < 50 pc sample, a 
six planet correction to the lowest mass bin (see Table 
12- footnote b). Fig. IHl shows a projection of Fig. \S\ 
onto the period axis. Planets are more clumped towards 
higher periods than are stellar companions. The Doppler 
planet detection method is not biased against short pe- 
riod planets. The Doppler stellar companion detections 
are not significantly biased for shorter periods or against 
longer periods in our samples analysis range (period < 5 
years) since Doppler instruments of much lower precision 
than those used to detect exoplanets are able to detect 
any Doppler companions of stellar mass. Thus this rep- 
resents a real difference in period distributions between 
stellar and planetary companions. 

The companions in Fig. [S] all have radial veloc- 
ity (Doppler) solutions. Some of the companions also 
have additional photometric, interferometric, astromet- 
ric or visual solutions. The exoplanet Doppler or- 
bits a re taken from the Extrasolar Planets Catalog 
((Schneider 2005). Only the planet orbiting the star 
HIP 108859 (HD 209458) has an additional photomet- 
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Fig. 5. — Brown Dwarf Desert in Mass and Period. Estimated 
companion mass M2 versus orbital period for the companions to 
Sun-like stars of our two samples: companions with hosts closer 
than 25 pc (large symbols) and those with hosts closer than 50 pc, 
excluding those closer than 25 pc (small symbols). The companions 
in the thick solid rectangle are defined by periods P < 5 years, and 
masses 10"^ < M2 ^ Mq, and form our less-biased sample of com- 
panions. The stellar (open circles), brown dwarf (grey circles) and 
planetary (filled circles) companions are separated by dashed lines 
at the hydrogen and deuterium burning onset masses of 80 Mj^p 
and 13 Mj^p respectively. This plot clearly shows the brown dwarf 
desert for the P < 5 year companions. Planets are more frequent at 
larger periods than at shorter periods (see Fig. The "Detected" , 
"Being Detected" and "Not Detected" regions of the mass-period 
space show the extent to which the high precisio n Doppler metho d 
is currently able to find companions fLincwcaver fc Gretherl200.'^) . 
see Appendix for discussion of M2 mass estimates. 

ric solution but this companion falls outside our less- 
biased region (M2 < Mj^p)- For the stellar com- 
panion data, the single- lined (SBl) and double-lined 
(SB2) spectroscopic binary orbits are primarily from 
the Ninth Catalogue of Spectroscopic Binary Orbits 
ijPourbaix et al\ 120041) with additional interferometric, 
astrometric or visual solutions from the 6th Catalog of 
Orbits of Visual Binary Stars (Washington Double Star 
Catalog, Hartkop f fc Mason 2004 ). Ma ny additional 
SB Is come from iHalbwachs et al\ l)2003(l . Stellar bi- 
naries a nd orbital solutio n s also c ome from lEndl et al.\ 
ll200l: IHalbwachs aO (l2000l): Maaeh. et al\ ll2003D: 



\T\xm&f et aV. ('2001'): 'Jones era/.l l|2002|) : iVogt a^J 
I2OO2); Zuckcr & Mazch (2001a). 

We examine the inclination distribution for the 30 
Doppler companions (d < 50 pc) with an astrometric or 
visual solution. We find that 24 of these 30 companions 
have a minimum mass larger than 80Mj„p (Doppler stel- 
lar candidates) and that 6 of these 30 companions have 
a minimum mass between ISMj^p and 80Mj„p (Doppler 
brown dwarf candidates) . These 6 Doppler brown dwarf 
candidates are a subset of the 16 Doppler brown dwarf 
candidates in the far sample that have an astrometric 
orbit derived with a confidence level greater than 95% 
from Hippa rcos measurements (Halbwa chs et al\ 120001: 
iZucker fcMazeh 2001a.) and are thus assumed to have 
an astrometric orbit. 
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Fig. 6. — Projection of Fig. Isl onto the period axis for the 25 
pc (dark grey) and 50 pc (light grey) samples. Planets are more 
clumped towards higher periods than are stellar companions. This 
would be a selection effect with no significance if the efficiency 
of finding short period stellar companions with the low precision 
Doppler technique used to find spectroscopic binaries, was much 
higher than the efficie ncy of finding e xopl anets with h igh precision 
spectroscopy. Konac ki et al\ 120041) and IPont et ali (2004) con- 
clude that the fact that the transit photometry method has found 
planets in sub 2.5 day periods (while the Doppler method has found 
none) is due to higher efficiency for small periods and many more 
target stars and thus that these two observations do not conflict. 
Thus there seems to be a real difference in the period distributions 
of stellar and planetary companions. 

As shown in Fig. [T] the inclination distribution is ap- 
proximately random for the 24 companions with a min- 
imum mass in the stellar regime whereas it is biased to- 
wards low inclinations for the 6 companions in the brown 
dwarf regime. All 6 of the Doppler brown dwarf candi- 
dates with an astrometric determination of their incli- 
nation have a true mass in the stellar regime. This in- 
cludes all 3 of the Doppler brown dwarf candidates that 
are companions to stars in our close sample (d < 25 pc) 
thus leaving an empty brown dwarf regime. Also shown 
in Fig. [71 is the distribution of the maximum values of 
sin{i) that would put the true masses of the remaining 
10 Doppler brown dwarf candidates with unknown incli- 
nations in the stellar regime. This distribution is sub- 
stantially less-biased than the observed sin{i) distribu- 
tion, strongly suggesting that the remaining 10 Doppler 
brown dwarf candidates will also have masses in the stel- 
lar regime. Thus astrometric corrections leave us with no 
solid candidates with masses in the brown dwarf regime 
from the 16 Doppler brown dwarf candidates in the far 
sample (d < 50 pc), consistent with the result obtained 
for the close sample. 

The size of the 25 pc and 50 pc samples, the extent 
to which they are being targeted for planets, and the 
number and types of companions found along with any 
associated corrections are summarised in Table |5J For 
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TABLE 2 

HippARCOS Sample, Doppler Targets and Detected Companions for Near and Far Samples 



Sample 


Hipparcos 


Doppler 


Companions 




Number 


Targ 


ct 


Total 


Planets 


BDs 


St 


ars 








Number 


%" 








Total 


SBl 


SB2 


d < 25 pc 


1509 


627 


42% 




22 










Sun-like 


464 


384 


83% 


59 (-f 15)^ 


19 ( + l\+4'=) 





40 (+7'', +3=) 


25 (9)'' 


15 (8)'' 


Dec < 0° 


211 


211 


100% 


20 (4-10)^ 


10 





10 (+7'^, -1-3°) 


8 (3)** 


2 (l)*" 


Dec > 0° 


253 


173 


68% 


39 


9 





30 


17 (6)'* 


13 (7)'' 


d < 50 pc 


6924 


2351 


34% 




58 










Sun-like 


3296 


2001 


61% 


198 (+80)-'' 


54 (^-6^-(-19=) 


19 


143 (-1-14'', -1-41") 


90 (18)*" 


53 (12)'* 


Dec < 0° 


1647 


1525 


93% 


72 (+74)^ 


33 ( + 19) = 





39 ( + 14"', -1-41") 


27 (7)'* 


12 (2)'' 


Doc > 0° 


1649 


476 


29% 


126 


21 


1 


104 


63 (11)'' 


41 (lO)*" 



" Percentage of Hipparcos stars that arc Doppler targets. 

Completeness correction in the lowest mass bin for the lower right corner of our sample in Fig, 5 lying in the "Being Detected" region (see 
Lineweaver, Grether Sz Hidas 2003), 

" Renormalization for planetary target population (384) being less than stellar companion mother population (464) (sec discussion in Section 2,1), 
^ Correction based on the most likely scenario that the southern stellar companions from Jones et al. 12002) have periods < 5 years, 
" Correction for north/south declination asymmetry in companion fraction after correcting for ijones et all I2003l detections (see Section 2.2). 
^ Total of corrections b through to c. 

9 Result from assuming < sin{i) >— 0.785 when i is unknown (see caption of Fig. 7 and Appendix). 

^ Number of these spectroscopic binaries with an additional astromctric or visual solution (sec Appendix). 



the stars closer than 25 pc, 59 have companions in the 
less-biased region (rectangle circumscribed by thick line) 
of Fig. Of these, 19 are exoplanets, are brown dwarfs 
and 40 are of stellar mass. Of the stellar companions, 25 
are SBls and 15 are SB2s. For the stars closer than 50 
pc, 198 have companions in the less-biased region. Of 
these, 54 are exoplanets, 1 is a brown dwarf and 143 are 
stars. Of the stellar companions, 90 are SBls and 53 are 
SB2s. 

We find an asymmetry in the north/south declination 
distribution of the Sun-like stars with companions, prob- 
ably due to undetected or unpublished stellar compan- 
ions in the south. The number of hosts closer than 25 
pc with planetary or brown dwarf companions are sym- 
metric in north/south declination to within one sigma 
Poisson error bars, but because more follow up work has 
been done in the north, more of the hosts with stellar 
companions with orbital solutions are in the northern 
hemisphere (30) compared with the southern (10). A 
comparison of our northern sample of hosts with stellar 
companions to the similarl y selected appr oximately com- 
plete sample of iHalbwachs et al\ l)2003(l indicates that 
our 25 pc northern sample of hosts with stellar compan- 
ions is also approximately complete. Under this assump- 
tion, the number of stellar companions missing from the 
south can be estimated by making a minimal correction 
up to the one sigma error level below the expected num- 
ber, based on the northern follow up results. Of the 464 
Sun-like stars closer than 25 pc, 211 have a southern dec- 
lination (Dec < 0°) and 253 have a northern declination 
(Dec > 0°) and thus - 25(25/211 « 30/253) stars in the 
south should have a stellar companion when fully cor- 
rected or 20 if we make a minimal correction. Thus we es- 
timate that we are missing at least ~ 10(= 20—10) stellar 
co mpanions in t he sou th, 7 of which have been detected 
by iJones et'of] l)2002() under the plausible assumption 
that the orbital periods of the companions detected by 
|jones et al. { 2002) are less than 5 yea rs. Alth ough these 
7 SBl stellar companions detected bv lJones et all lf200^ 
have as yet no published orbital solutions, we assume 
that t he SBl stellar companions detected bv lJones et all 
l|2002D have P < 5 years since they have been observed as 



part of the high Doppler precision program at the Anglo- 
Australian Observatory (started in 1998) for a duration 
of less than 5 years before being announced. The addi- 
tional estimated stellar companions are assumed to have 
the same mass distribution as the other stellar compan- 
ions. 

We can similarly correct the declination asymmetry in 
the sample of Sun-like stars closer than 50 pc. We find 
that there should be, after a minimal correction, an ad- 
ditional 55 stars that are stellar companion hosts in the 
southern hcmipshere. 14 of these 55 stellar companion s 
are assumed to have been detected bv lJones' 
An asymmetry found in the planetary companion frac- 
tion in the 50 pc sample due to the much larger number 
of stars being monitored less intensively for exoplanets 
in the south (~ 2% = 33/1525) compared to the north 
(~ 4% — 21/476) results in a correction of 19 planetary 
companions in the south. The results given in Table|21are 
done both with and without the asymmetry corrections. 

Unlike the 25 pc sample for which we are confident 
that the small corrections made to the number of com- 
panions will result in a reliable estimate of a census, cor- 
recting the 50 pc sample for the large number of missing 
companions is less reliable. This is so because if it were 
complete, the 50 pc sample would have approximately 
8 times the number of companions as the 25 pc sample, 
since the 50 pc sample has 8 times the volume of the 25 pc 
sample. However, the incomplete 50 pc sample has only 
~ 7{— 3296/464) times the number of Hipparcos stars, 
~ 5(= 2001/384) times as many exoplanet targets and 
~ 3 times as many companions as the 25 pc sample. Thus 
rather than correcting both planetary and stellar com- 
panions by large amounts we show in Section 3 that the 
relative number and distribution of the observed plane- 
tary and stellar companions (plus a small completeness 
correction for the "Being Detected" region of 6 planets 
and an additional 14 probable stellar companions from 
[^nes et al. (20Q2) - see Table |21l remains approximately 
unchanged when compared to the corrected companion 
distribution of the 25 pc sample. Analyses both with 
and without a correction for the north/south asymmetry 
produce similar results for the brown dwarf desert (Table 
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Fig. 7. — Astrometric inclination distribution for close com- 
panions {d < 50 pc) with a minimum mass larger than 80Mj„p 
(Doppler stellar candidates - TOP) and between ISMj^p and 
80Mj„p (Doppler brown dwarf candidates - BOTTOM). There are 
24 companions with astrometric solutions and a minimum mass in 
the stellar regime and 6 with a minimum mass in the brown dwarf 
regime. The inclination distribution is approximately random for 
companions with a minimum mass in the stellar regime whereas it is 
biased towards low inclinations for companions in the brown dwarf 
regime. All 6 astrometric determinations of sin{i) for brown dwarf 
candidates put their true mass in the stellar regime. Also shown is 
the distribution of the maximum values of sin(i) that would place 
the true masses of the remaining 10 brown dwarf candidates with- 
out astrometric or visual solutions in the stellar regime. A distri- 
bution less-biased than the observed sinii) distribution would be 
required. This strongly suggests that the 10 candidates without 
astrometric or visual solutions will also have masses in the stel- 
lar regime. Therefore, astrometric corrections leave us with no 
solid candidates with masses in the brown dwarf region. Two weak 
brown dwarf candidates are worth mentioning. HD 114762 has 
a minimum mass below 13Mj„p. However, to convert minimum 
mass to mass, we have assumed random inclinations and have used 
< sinii) >Kd 0.785. This conversion puts the estimated mass of HD 
114762 in the brown dwarf regime (M Z 13Mj„p). In Fig. El this 
is the only companion lying in the brown dwarf regime. Another 
weak brown dwarf candidate is the only candidate that requires a 
sin{i) < 0.2 to place its mass in the stellar regime. 

13. 

3. COMPANION MASS FUNCTION 

The close companion mass function to Sun-like stars 
clearly shows a brown dwarf desert for both the 25 pc 
(Fig. |HJ and the 50 pc (Fig. |5J samples. The numbers 
of both the planetary and stellar mass companions de- 
crease toward the brown dwarf mass range. Both plots 
contain the detected Doppler companions, shown as the 
grey histogram, within our less-biased sample of compan- 
ions (P < 5 years and M2 > 10~^Mq, see Section 2.2). 
The hatched histograms at large mass show the subset of 
the stellar companions that are not included in any of the 
exoplanet Doppler surveys. A large bias against stellar 
companions would have been present if we had only in- 
cluded companions found by the exoplanet surveys. For 
multiple companion systems, we select the most massive 
companion in our less biased sample to represent the sys- 
tem. We put the few companions (3 in the 25 pc sample, 
6 in the 50 pc sample) that have a mass slightly larger 



than 1 Mq in the largest mass bin in the companion mass 
distributions. 

Fitting straight lines using a weighted least squares 
method to the 3 bins on the left-hand side (LHS) and 
right-hand side (RHS) of the brown dwarf region of the 
mass histograms (Figs. |H1 & El, gives us gradients of 
-15.2 ± 5.6 (LHS) and 22.0 ± 8.8 (RHS) for the 25 pc 
sample and -9.1 ± 2.9 (LHS) and 24.1 ± 4.7 (RHS) for 
the 50 pc sample. Since the slopes have opposite signs, 
they form a valley which is the brown dwarf desert. The 
presence of a valley between the negative and positive 
sloped lines is significant at more than the 3 sigma level. 
The ratio of the corrected number of companions in the 
less-biased sample on the LHS to the RHS along with 
their poisson error bars is (24±9)/(50±13) = 0.48±0.22 
with no companions in the middle 2 bins for the 25 pc 
sample. For the larger 50 pc sample the corrected less- 
biased LHS/RHS ratio is (60 ± 14)/(157 ± 22) = 0.38 ± 
0.10, with 1 brown dwarf companion in the middle 2 bins. 
Thus the LHS and RHS slopes agree to within about 1 
sigma and so do the LHS/RHS ratios, indicating that the 
companion mass distribution for the larger 50 pc sample 
is not significantly different from the more complete 25 pc 
sample and that the relative fraction of planetary, brown 
dwarf and stellar companions is approximately the same. 
A comparison of the relative number of companions in 
each bin in Fig. IHlwith its corresponding bin in Fig. El 
produces a best-fit of = 1-9. 

To find the driest part of the desert, we fit separate 
straight lines to the 3 bins on either side of the brown 
dwarf desert (solid lines) in Figs. |S1 & El The deepest 
part of the valley where the straight lines cross beneath 
the abscissa is at M = SitllMjup and M = 43 taa^Jup 
for the 25 and 50 pc samples respectively. These results 
are summarized in Table O The driest part of the desert 
is virtually the same for both samples even though we 
see a bias in the stellar binarity fraction of the 50 pc 
sample (Fig. EJ. We have done the analysis with and 
without the minimal declination asymmetry correction. 
The position of the brown dwarf minimum and the slopes 
are robust to this correction (see Table EJ. 

The smaller 25 pc Sun-like sample contains 464 stars 
with 16.0%±5.2% of these having companions in our cor- 
rected less-biased sample. Of these ~ 16% with compan- 
ions, 5.2%±1.9% are of planetary mass and 10.8%±2.9% 
are of stellar mass. None is of brown dwarf mass. This 
agrees with previ ous estimates of stellar binarity such 
as that found bv iHalbwachs et al] (|2003D of 14% for a 
sample of G-dwarf companions with a slightly larger pe- 
riod range (P < 10 years). The planet fr action agrees 
with t he fraction 4%±1% found in lLineweav er & Grcthe3 
()2003ft when most of the known exoplanets are consid- 
ered. The 50 pc sample has a large incompleteness due 
to the lower fraction of monitored stars (Fig. but as 
shown above, the relative number of companion planets, 
brown dwarfs and stars is approximately the same as for 
the 25 pc sample. The 50 pc sample has a total com- 
panion fraction of 15.6% ± 2.8%, where 4.3% ± 1.0% of 
the companions are of planetary mass, are of 

brown dwarf mass and 11.2% ± 1.6% are of stellar mass. 
Table 01 summarizes these companion fractions. 

Surveys of the multiplicity of nearby Sun-like stars 
yield the relative numbers of single, double and multiple 
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Companion Moss M2 [M©] 

Fig. 8. — Brown Dwarf Desert in Close Sample. Histogram of 
the companions to Sun-like stars closer than 25 pc plotted against 
mass. The grey histogram is made up of Doppler detected compan- 
ions in our less-biased (P < 5 years and M2 > IO^'^Mq) sample. 
The corrected version of t his less-biased sam ple includes an extra 7 
probable SBl stars from IJones et a;j|20 ^ (Tablel^- footnote d) 
and an extra 3 stars from an asymmetry in the host declination dis- 
tribution (Table 121- footnote e). The planetary mass companions 
are also renormalized to account for the small number of Hipparcos 
Sun-like stars that are not being Doppler monitored (21% renor- 
malization, Table 1^- footnote c) and a 1 planet correction for the 
undersampling of the lowest mass bin due to the overlap with the 
"Being Detected" region (Tablel^l- footnote b). The hatched his- 
togram is the subset of detected companions to hosts that are not 
included on any of the exoplanet search target lists and hence shows 
the extent to which the exoplanet target lists are biased against the 
detection of stellar companions. Since instruments with a radial 
velocity sensitivity Kg < 40 m/s (see Eq. 2 of Appendix) were 
used for all the companions, we expect no other substantial bi- 
ases to affect the relative amplitudes of the stellar companions on 
the right-hand side (RHS) and the planetary companions on the 
left-hand side (LHS). The brown dwarf mass range is empty. 



star systems. According to lDuauennov &: Mavorl l) 19911) . 

51% of star systems are single stars, 40% are double star 
systems, 7% are triple and 2% are quadruple or more. 
Of the 49%(= 40-1-7-1-2) which are stellar binaries 
or multiple star systems, 11% have stellar companions 
with periods less than 5 years and thus we can infer that 
the remaining 38% have stellar companions with P > 5 
years. Among the 51% without stellar companions, we 
find that ^ 5% have close (P < 5 years) planetary com- 
panions with 1 < A'l/Mjup < 13, while < 1% have close 
brown dwarfs companions. 

The Doppler method should preferentially find planets 
around lower mass stars where a greater radial velocity 
is induced. This is the opposite of what is observed as 
shown in Figs. 1101 and ITTl where we split the 25 and 50 
pc samples respectively into companions to hosts with 
masses above and below 1 Mq. We scale these smaller 
samples to the size of the full 25 and 50 pc samples 
(Figs. IHl and El respectively). The Doppler tech nique 
is also a function of P — color ijSaar et oi.lll998j) with 
the level of systematic errors in the radial velocity mea- 
surements, decreasing as we move from high mass to low 
mass {B — V — 0.5 to P — ^ = 1.0) through our two 
samples, peaking for late K spectral type stars before in- 
creasing for the lowest mass M type stars again. Hence 




10-2 10- 

Componlon Moss M2 [M©] 

Fig. 9. — Same as Fig. ISlbut for the larger 50 pc sample renor- 
malized to the size of the 25 pc sample. Fitting straight lines using 
a weighted least squares fit to the 3 bins on the LHS and RHS, 
gives us gradients of —9.1 ± 2.9 and 24.1 ± 4.7 respectively (solid 
lines). Hence the brown dwarf desert is significant at more than the 
3 sigma level. These LHS and RHS slopes agree to within about 1 
sigma of those in Fig. O The ratio of the number of companions 
on the LHS to the RHS is also about the same for both samples. 
Hence the relative number and distribution of companions is ap- 
proximately the same as in Fig. O The separate straight line fits 
to the 3 bins on the LHS and RHS intersect at M = 43 'tl^t-^^J^p 
beneath the abscissa. Approximately 16% of the stars have com- 
panions in our less-biased region. Of these, 4.3% ± 1.0% have com- 
panions of planetary mass, 0.1^q'^% have brown dwarf companions 
and 11.2%±1.6% have companions of stellar mass. We renormalize 
the mass distribution in this figure by comparing each bin in this 
figure with its corresponding bin in Fig. jSland scaling the vertical 
axis of Fig. j^lso that the difference in height between the bins is 
on average a minimum. We find that the optimum renormalization 
factor is 0.33. This plot does not include the asymmetry correction 
for the planetary and stellar companions discussed in Section 2.2 
and shown in Table 1^ 

again finding planets around the lower mass stars (early 
K spectral type) in our sample should be easier. 

3.1. Comparison with Other Results 

Although there are s ome similariti e s, the companion 
mass function found bv 'HcacroJ l)1999D : IZucker fc Mazehl 
(2001b); Mazeh et al. (2003) is different from that shown 
in Figs. IHl & IHl Our approach was to normalize the com- 
panion numbers to a well-defined sub-sample of Hippar- 
cos stars whereas these authors use two different sam- 
ples of stars, one to find the planetary companion mass 
function and another to find the stellar companion mass 
function, which are then normalized to each other. The 
different host star properties and levels of completeness 
of the two samples may make this method more prone 
than our method, to biases in the frequencies of compan- 
ions. 

l i. r ■ iHeacoxl l(T999|) and iZucker fc Ma^ l)2001bl) 
combined the companions o f the stellar mass sample of 
iDuauennov fc MavoH l|199lD with the known substellar 
companions, but identified different mass functions for 
the planetary mass regime below 10 Mj^p and similar 
flat distributions in logarithmic mass for b rown dwarf 
and stellar mass companions. „,Heaco?d l|1999D found that 
the logarithmic mass function in the planetary regime is 
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TABLE 3 

Companion Slopes and Companion Desert Mass Minima 



Sample 




Asymmetry 
Correetion 


Figure 


LHS slope 


RHS slope 


Slope Minima^ 

[A/.7„p] 


d < 25 pc 




Yes 


8 


-15.2 ± 5.6 


22.0 ± 8.8 


+5S 
■^^ -18 


d < 25 pc 




No 




-15.2 ± 5.6 


20.7 ± 8.5 


30 


d < 50 pc 




Yes 




-9.4 ± 3.0 


24.3 ± 4.6 




d < 50 pc 




No 


9 


-9.1 ± 2.9 


24.1 ± 4.7 


43 1^3 


d < 25 pc ,S 


i Ml < IMq 


Yes 


10 


-17.5 ± 5.4 


19.4 ± 10.7 


18|^f 


d < 50 pc ^ 


i Ml < 1M0 


No 


11 


-5.9 ± 5.1 


25.2 ± 11.4 


39 —23 


d < 25 pc ,S 


i Ml > IMq 


Yes 


10 


-12.4 ± 9.2 


20.0 ± 10.9 


-26 


d < 50 pc .S 


L Ml > lAfg 


No 


11 


-12.2 ± 8.2 


21.1 ± 10.4 


45 

-21 



" values of mass where the best-fitting hnes, to the LHS and RHS, intersect. The errors given are from the range between the two intersections 
with the abscissa. 



TABLE 4 
Companion Fraction Comparison 



Sample 


Asymmetry 
Correetion 


Figure 


Total % 


Planetary % 


Brown Dwarf % 


Stellar % 


d < 25 pe 
d < 25 pe 


Yes 
No 


8 


16.0 ± 5.2 
15.3 ± 5.0 


5.2 ± 1.9 
5.2 ± 1.9 


Q 0+0.4 
^■^ -0.0 

"■^ -0.0 


10.8 ± 2.9 
10.1 ± 2.7 


d < 50 pc 
d < 50 pc 


Yes 
No 


9 


15.6 ± 2.8 
15.6 ± 2.8 


4.4 ± 1.0 
4.3 ± 1.0 


1+"-^ 

^'■^ -0.1 

1 +"-2 

^■^ -0.1 


11.1 ± 1.6 

11.2 ± 1.6 


d < 25 pe & Ml < IMq 
d < 50 pe & Ml < IMq 


Yes 
No 


10 
11 


16.0 ± 5.8 
15.6 ± 6.0 


4.2 ± 1.9 
2.6 ± 1.7 


0+"'' 

^'■^ -0.0 

2+01 

^'■-^ -0.2 


11.8 ± 3.5 
12.8 ± 3.9 


d < 25 pe & Ml > IMq 
d < 50 pe & Afi > IMq 


Yes 
No 


10 
11 


16.0 ± 7.0 
15.6 ± 6.7 


6.6 ± 3.1 
6.2 ± 2.9 


0+"'' 

"■'J -0.0 

0+"'' 

^'■^ -0.0 


9.4 ± 3.5 
9.4 ± 3.4 
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Fig. 10. — Same as Fig. ISl but for the 25 pc sample split into 
companions to lower mass hosts (A/i < IA/0) and companions to 
higher mass hosts {M\ > IMq). The lower mass hosts have 4.2% 
planetary, 0.0% brown dwarf and 11.8% stellar companions. The 
higher mass hosts have 6.6% planetary, 0.0% brown dwarf and 
9.4% stellar companions. The Doppler method should preferen- 
tially find planets around lower mass stars where a greater radial 
velocity is induced. This is the opposite of what we observe. To 
aid comparison, both samples are scaled such that they contain the 
same number of companions as the full corrected less-biased 25 pc 
sample of Fig. |S] 

best fit by a power-law {d N / dlogM oc M^) with in dex F 
between and -1 wliereas |Zucker fc Mazebl (|2001b) find 
an approximately flat dist ribution (power-law with index 
0). Our work here and in iLineweaver fc Cxretheil lEffll 




10-2 10- 
Companlon Mass M2 [M©] 

Fig. 11. — Same as Fig. |5]but for the 50 pc sample split into 
companions to lower mass hosts (Mi < IMq) and companions to 
higher mass hosts (Afi > IMq). Both samples are scaled such that 
they contain the same number of companions as the corrected less- 
biased 50 pc sample of Fig. |^ Also shown are the linear best-fits 
to the planetary and stellar companions of the two populations. 



suggests that neither the stellar nor the planetary com- 
panion distributions are flat (F — —0.7). Rather, they 
both slope down towards the br own dwarf des ert, more 
in agreement with the results of lHeaco^ 1)19991) . 

The work mos t similar to ours is probably 
ijMazeh et al\ |2003) who looked at a sample of main 
sequence stars with primaries in the range 0.6 — 0.85 Mq 
and P < 3000 days using infrared spectroscopy and 
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Fig. 12. — The mass function of companions to Sun-like 
stars (lower left) compared to the initial mass function (IMF) 
of cluster stars (upper right). Our mass function of the com- 
panions to Sun-like stars is shown by the green dots (bigger 
dots are the d < 25 pc sample, smaller dots are the d < 50 
pc sample). The linear slopes we fit to the data in Fig. |S] 
are also shown along with their error. Data for the number 
of stars and brown dwarfs in the Orion Nebula Cluster (ONC) 
(circles), Pleiades cluster (triangles) and M35 cluster (squares ) 
come fi-om Hillenbrand & Carocnter (2000); Slcs nick ef oD 120(111 . 
IMoraux et al. 1 2003) and Barrado v Navascu cs^^an ~ l20Ml) re- 
spectively and are normalized such that they overlap for masses 
larger than IMq where a single power-law slope applies. The ab- 
solute normalization of cluster stars is arbitrary, while the com- 
panion mass function is normalized to the IMF of the cluster stars 
by scaling the three companion points of stellar mass to be on av- 
er age ~ 7% for P < 5 y ears (derived from the stellar multiplicity 
of IDuguennov &: Mavon Il99ll) discussed in Section 3, combined 
with our estimate that 11% of Sun-like stars have stellar secon- 
daries). The average power-law IMF derived from various values 
of the slope of the IMF quoted in the literature ( Hillenbrand 2003) 
is shown as larger red dots along with two thin red lines showing 
the root-mean-square error. If the turn down in the number of 
brown dwarfs of the IMF is due to a selection effect because it is 
hard to detect brown dwarfs, then the two distributions are even 
more different from each other. For clarity the smaller green dots 
are shifted slightly to the right. 

combined them with the known substellar companions of 
these main sequence stars and found that in logarithmic 
mass the stellar companions reduce in number towards 
the brown dwarf mass range. This agrees with our 
results for the shape of the stellar mass companion 
function. However, they identify a flat distribution 
for the planetary mass companions in contrast to our 
non-zero slope (see Table EJ. Mazeh et al. ( 2003) found 
the frequency of stellar and planetary companions 
(M2 > 1 Mjup) to be 15% (for stars below 0.7 Mq) and 
3% respectively. This compares with our estimates of 
8% (for stars below 0.7 Mr.-)) and 5%. The larger period 
range used by iMazeh e t al. ( 2003) can account for the 
difference in stellar companion fractions. 

4. COMPARISON WITH THE INITIAL MASS FUNCTION 

Brown dwarfs found as free-floating objects in the so- 
lar neighbourhood and as members of young star clus- 
ters have been used to extend the initial mass function 
(IMF) well into the brown dwarf regime. Comparing the 
mass function of our sample of close-orbiting compan- 
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Fig. 13. — The initial mass function (IMF) for clusters repre- 
sented by a series of power-law slopes (Hillenbrand 2003). Each 
point represents the power-law slope claimed to apply within the 
mass range indicated by the horizontal lines. Although the IMF 
is represented by a series of power-laws, the IMF is not a power- 
law for masses less than IMq where the slope continually changes. 
The green dots show the slope of the companion mass function to 
Sun-like stars between the bins of Figs. |S]&|5]with the larger and 
smaller dots respectively. The linear fits to the data in Fig. Island 
their associated error are shown by the curves inside the grey re- 
gions. The power-law fit of Lineweaver & Grether (2003) (shown 
as the green dot with a horizontal line indicating the range over 
which the slope applies) is consistent with these fits. The larger 
red dots with error bars represent the average power-law IMF with 
a root-mean-square error. F and —a are the respective logarith- 
mic and linear slopes of the mass function. The logarithmic mass 
power-law distribution is dN/dlogM oc and the linear mass 
power-law distribution is dN/dM oc M~" where F = 1 — a. The 
errors on the fits of Fig. |S] get smaller at M ~ 10^'^ Mq and 
Af ~ 1 Mq since as log(M/MQ ) tends to ±00, F tends to 0. This 
can also be seen in Fig. I12l where the slopes of the upper and lower 
contours become increasingly similar. 

ions of Sun-like stars to the IMF of single stars indicates 
how the environment of a host affects stellar and brown 
dwarf formation and/or migration. Here we quantify how 
different the com panion mass function is from the IMF 
ijHalbwachs et aLli200(]() . 

The galactic IMF appears to be remarkably universal 
and independent of environment and metallicity with the 
possible exception of the substellar mass regime. A weak 
empirical trend with metallicity is suggested for very low 
mass stars and brown dwarfs where more metal rich en- 
vironme nts mav be pr oducing relatively more low mass 
objects llKrouDall2002|) . This is consistent with an ex- 
trapolation up in mass from the trend found in exoplanet 
hosts. The IMF is often represented as a power-law, al- 
though this only appea rs to be accurate for stars with 
masses above ~ IMq l|Hillenbrandl 12003(1 . The stellar 
IMF slope gets flatter towards lower masses and extends 
smoothly and continously into the substellar mass regime 
where it appears to turn over. 

Free floating brown dwarfs may be formed either as 
ejected stellar embryos or from low mass protostellar 
cores that have lost their accretion envelopes due to 
photo-evaporation from the chance proximity of a nearby 
massive star (Kroupa & Bouvier 2003). This hypothe- 
sis may explain their occurence in relatively rich star 
clusters such as the Orion Nebula cluster and their vir- 
tual absence in pre-main sequence stellar groups such as 
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Taurus- Auriga. 

In Figs. El & El we compare the mass function of 
companions to Sun-like stars with the IMF of cluster 
stars. The mass function for companions to Sun-like stars 
is shown by the green dots from Figs. IHlandlHl (bigger dots 
are the c? < 25 pc sample and smaller dots are the d < 50 
pc sample). The linear slopes from Fig. |S1 and their 
one sigma confidence region are also shown. Between 
log{M/MQ) K. -1.0 and -0.5 (O.IMq < M < 0.3Mq) 
the slopes are similar. However, above O.SM© and below 
O.IMq the slopes become inconsistent. Above O.SM© the 
slopes, while of similar magnitude are of opposite sign 
and below O.lM© the companion slope is much steeper 
than the IMF slope. The IMF for young clusters (yellow 
dots) is statistically indistinguishable from that of older 
stars (blue dots) and follows the average IMF. 

5. SUMMARY AND DISCUSSION 

We analyse the close-orbitting (P < 5 years) plane- 
tary, brown dwarf and stellar companions to Sun-like 
stars to help constrain their formation and migration sce- 
narios. We use the same sample to extract the relative 
numbers of planetary, brown dwarf and stellar compan- 
ions and verify the existence of a brown dwarf desert. 
Both planetary and stellar companions reduce in num- 
ber towards the brown dwarf mass range. We fit the 
companion mass function over the range that we anal- 
yse (0.001 < M/Mq < 1.0) by two separate straight fines 
fit separately to the planetary and stellar data points. 
The straight lines intersect in the brown dwarf regime, 
at M = 31 Mjup- This result is robust to the decli- 
nation asymmetry correction (Table |2J). 

The period distribution of close-orbitting (P < 5 years) 
companion stars is different from that of the planetary 
companions. The close-in stellar companions are fairly 
evenly distributed over logP with planets tending to be 
clumped towards higher periods. We compare the com- 
panion mass function to the IMF for bodies in the brown 
dwarf and stellar regime. We find that starting at 1 Mq 
and decreasing in mass, stellar companions continue to 
reduce in number into the brown dwarf regime, while 
cluster stars increase in number before reaching a maxi- 
mum just before the brown dwarf regime (Fig. 13). This 
leads to a difference of at least 1.5 orders of magnitude 
between the much larger number of brown dwarfs found 
in clusters to those found as close-orbitting companions 
to Sun-like stars. 

The period distribution of close-orbiting companions 
may be more a result of post-formation migration and 
gravitational jostling than representive of the relative 
number of companions that are formed at a specific dis- 
tance from their hosts. The companion mass distribution 
is more fundamental than the period distribution and 
should provide better constraints on formation models, 
but our ability to sample the mass distribution is only 
for P < 5 years. 

We show in Figs. ^1 and that lower mass hosts 
have more stellar companions and fewer giant planet 
companions while higher mass hosts have fewer stellar 
companions but more giant planet companions. The 
brown dwarf desert is generally thought to exist at close 
separations ^ 3 AU (or equivalently P < 5 years) 
llMarcv fc Butleill20qol) but may disappear at wider sep- 
arations. iGizis et all l|200ij) suggests that at very large 



separations (> 1000 AU) brow n dwarf companions may 
be mo re common. However, iMcCarthv fc ZuckermanI 
(j2004|) in their observation of 280 GKM stars find only 1 
brown dwarf between 75 and 1200 AU. Gi zis et al\ l|2003D 
reports that 15% ± 5% of M/L dwarfs are brown dwarf 
binaries with separations in the range 1.6 — 16 AU. This 
falls to 5% ± 3% of M/L dwarfs with separations less 
than 1.6 AU and none with separations greater than 16 
AU. This differs greatly from the brown dwarfs orbiting 
Sun- like stars but is consistent with our host /niinimum- 
companion-mass relationship, i.e., we expect no short pe- 
riod brown dwarf desert around M or L type stars. 

Three systems containing both a companion with a 
minimum mass in the planetary regime and a com- 
panion with a minimum mas s in the bro wn dwarf 
regime are known - HP 168443 llMarcv eit ftj! jl20011 . HD 
202206 iftldTv oi.|l2?)0llCorreia et QiJl2004D and GJ 86 
(jOueloz et a/.ll2000t lEls et a^.ll200H) ~Our analysis sug- 
gests that both the Msm(i)-brown dwarfs orbiting HD 
168443 and HD 202206 are probably stars (see Section 
2.2 for our false positive brown dwarf correction). If the 
Msm(i)-planetary companions in these 2 systems are 
coplanar with the larger companions then these "plan- 
ets" may be brown dwarfs or even stars. GJ 86 con- 
tains a possible brown dwarf detected orbiting at ~ 20 
AU (P > 5 years) and so was not part of our analysis. 
However this does suggest that systems containing stars, 
brown dwarfs and planets may be possible. 

We find that approximately 16% of Sun-like stars have 
a close companion more massive than Jupiter. Of these 
16%, 11% ± 3% are stellar, < 1% are brown dwarf and 
5% lb 2% are planetary comp anions (Table Although 
iLineweaver fc Gretheil l|2003f l show that the fraction of 
Sun-like stars with planets is greater than 25%, this is for 
target stars that have been monitored the longest (~ 15 
years) and at optimum conditions (stars with low-level 
chromospheric activity or slow rotation) using the high 
precision Doppler method. W hen we limit the analysis of 
ILineweaver fc GretheJ ll200l to planetary companions 
with periods of less than 5 years and masses larger than 
Jupiter, we find the same value that we calculate here. 
When we split our sample of companions into those with 
hosts above and below IMq, we find that for the lower 
mass hosts: 11.8% have stellar, < 1% have brown dwarf 
and 4.2% have planetary companions and that for the 
higher mass hosts: 9.4% have stellar, < 1% have brown 
dwarf and 6.6% have planetary companions respectively 
(Table . More massive hosts have more planets and 
fewer stellar companions than less massive hosts. These 
are marginal results but are seen in both the 25 and 50 
pc samples. 

The constraints that we have identified for the com- 
panions to Sun-like stars indicate that close orbitting 
brown dwarfs are very rare. The fact that there is a 
close-orbitting brown dwarf desert but no free fioating 
brown dwarf desert suggests that post-collapse migration 
mechanisms may be responsible for this relative dearth 
of observable brown dwarfs rather than some intrinsic 
minimum in fragmentation an d gravitational co llapse in 
the brown dwarf mass regime ilda fc Lin 112001 . What- 
ever migration mechanism is responsible for putting hot 
Jupiters in close orbits, its effectiveness may depend on 
the mass ratio of the object to the disk mass. Since there 
is evidence that disk mass is correlated to host mass, the 
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as 



migratory m echanism may be correlated to host mass, 
proposed bv lArmitage fc Bonnelll l|2002D . 
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7. APPENDIX: COMPANION MASS ESTIMATES 

The Doppler method for companion detection cannot 
give us the mass of a companion without some additional 
astrometric or visual solution for the system or by mak- 
ing certain assumptions about the unknown inclination 
except in the case where a host star and its stellar com- 
panion have approximately equal masses and a double- 
lined solution is available. Thus to find the companion 
mass M2 that induces a radial velocity Ki in a host star 
of mass All we use (see Heacox 1999) 



,2^G,i/3 M2sin{i) 
^1 = i—fr) 



1 



P ' (Mi+M2)2/3(l_g2)l/2 



(1) 



This equation can be expressed in terms of the mass 
function /(m) 

/M = /.^ , .^^2 = ^ (2) 



(Ml + M^f 



2ttG 



Eq. 3 can then be expressed in terms of a cubic equa- 
tion in the mass ratio q = M2/M1, where Y = f{m)/Mi. 



For planets (Mi >> M2) we can simplify Eq. 2 and 
directly solve for M2sin{i) but this is not true for larger 
mass companions such as brown dwarfs and stars. We use 
[Cox (2000) to relate host mass to spectral type. When 
a double-lined solution is available, the companion mass 
can be found from q = M2/M1 = Ki/K2- 

For all single- lined Doppler solutions, where the incli- 
nation i of a companion's orbit is unknown (no astromet- 
ric or visual solution), we assume a random distribution 
P(i) for the orientation of the inclination with respect to 
our line of sight, 



P{i)di ~ sin{i)di 



(4) 



q^sin^{i) -Yq^ - 2Yq-Y = 



(3) 



From this we c an find proba bility distributions for 
sin{i) and sin^{i). iHeacoxl l)199,'Tj) and others suggest us- 
ing either the Richardson-Lucy or Mazeh- Goldberg algo- 
rith ms to approx imate the inclination distribution. How- 
ever, j^^gS^^] (j[r991() and lTrimbld l|199flll argue that for 
low number statistics, the simple mean method produces 
similar results to the more complicated methods. We 
have large bin sizes and small number statistics, hence 
we use this method. The average values of the sin(i) and 
sin^{i) distributions assuming a random inclination are 
< sin{i) >= 0.785 and < sin^{i) >= 0.589, which are 
used to estimate the mass for planets and other larger 
single-lined spectroscopic binaries respectively. For ex- 
ample, in Fig. O of the 198 mass estimates in the 
50 pc sample, 53 (27%) come from visual double-lined 
Doppler solutions, 6 (3% ) come from infrared double- 
hned Doppler solutions ijMazeh et all 12003(1 . 18 (9%) 
come from knowing the inclination (astrometric or visual 
solution also available for system), 10 (5%) come from as- 
suming that Doppler brown dwarf candidates have low 
inclinations, 55 (28%) come from assuming < sin(i) >= 
0.785 and 56 (28%) from assuming < sin^{i) >= 0.589. 
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